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Hydrodynamic Equations in Quantum Hall Systems at Large CurrentsQ 
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Hydrodynamic equations (HDEQs) are derived which describe spatio-temporal evolutions of 
the electron temperature and the chemical potential of two-dimensional systems in strong mag- 
netic fields in states with large diagonal resistivity appearing at the breakdown of the quantum 
Hall effect. The derivation is based on microscopic electronic processes consisting of drift motions 
in a slowly-fluctuating potential and scattering processes due to electron-electron and electron- 
phonon interactions. In contrast with the usual HDEQs, one of the derived HDEQs has a term 
with an energy flux perpendicular to the electric fleld due to the drift motions in the magnetic 
field. As an illustration, the current distribution is calculated using the derived HDEQs. 
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§1. Introduction 

The diagonal resistwity pxx is extremely small in the 
quantum Hall effectElcP (QHE) which occurs when the 
current density j is small. When j is increased up to a 
critical value, on the other hand, pxx increases by sev- 
eral orders of magnitude within a narrow range of j, 
and the QHE breaks downaQff (reviews are given in 
refs. Nm^hh Arrnmpanying the breakdown of the QHE, 
spatiaflBBBB) and temporaBBB evolutions of 
Pxx have been observed, and other various phenomena 
are expected such as pattern formations and oscilla- 
tory instabilities which piiave been extensively studied 
in bulk semiconductors. B* Unique phenomena are also 
expected in the present systems with strong magnetic 
fields, piYiuch rhfli^ been partly revealed in the experi- 
ments.BBy.BB) 

Theioipetically, on the other hand, a hot-electron 
modelBBB) has been proposed as a mechanism of the 
breakdown of the QHE. It is proposed in this model that 
the electron heating is responsible for the increase of pxx 
at the breakdown, that is, the electron temperature 
is the key variable in determining pxx- Later the hot- 
electron model has been supporleii hv^served spatial 
variations of pxx at a constant j Jj-tJ'BB which are un- 
derstood if we assume that spatial variations of are 
responsible for those of pxx ■ The importance of Te has al- 
ready been suggested by the previous experiments show- 
ing a bistability oi^^m&W-pxx state and a large- pa;^; state 
at the breakdown,!' dEP since the presence of such bista- 
bility requires another variable, in addition to j, to de- 
termine Pxx- Based on the hot-electron model, a hydro- 
dynamic equation (HDEQ) to calculate spatio-temporal 
variations of To has been proposed for one-dimensional 
variations in the previous work.B However, its valid- 
ity in the present quantum Hall systems is not known, 
because it was not derived on the basis of microscopic 

* to be published in J. Phys. Soc. Jpn. 71 (2002) No.l. 
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electronic processes. The theoretical framework to pre- 
dict a variety of phenomena in the present systems has 
not been established. 

In this study, we derive a set of HDEQs to describe 
spatio-temporal variations of the electron temperature 
Te and the chemical potential p in large-pa;a; states on 
the basis of microscopic processes in a two-dimensional 
electron system (2DES) in strong magnetic fields in the 
presence of a slowly-fluctuating potential due to ionized 
donors. As an example of its applications, we calculate 
the current distribution. One of the HDEQs in this paper 
is a generalization of the HDEQ given in our previous 
paper ,0 which is the equation for in the case of one- 
dimensional spatial modulations restricted to the current 
direction. 

The organization of this paper is as follows. In §2, 
three major assumptions employed in deriving the HD- 
EQs are introduced with several minor assumptions, and 
a model for electronic states and prpcesses, which has 
been proposed in our previous paper J2j) is described. In 
§3, two HDEQs are derived which describe evolutions of 
the electron density and the energy density, respectively. 
In §4, formulas of the electron and energy flux densities 
as well as the energy loss, appearing in the HDEQs, are 
derived. In §5, values of coefficients in the HDEQs are es- 
timated, and scaled HDEQs are introduced. The critical 
electric field is calculated using the estimated values of 
coefficients in steady and uniform states at even-integer 
filling factors, and is compared with observed values. In 
§6, the HDEQs are applied to obtain the current distri- 
bution at even-integer filling factors. In §7, conclusions 
and discussions are given. 

§2. Assumptions and Model 

2. 1 Single- Electron Approximation and Scattering Pro- 
cesses 

The potential Vdonor produced by ionized donors fiuc- 
tuates_in the plane of the 2DES. Its typical length 
scale,B or the average distance between a potential hill 
and a neighboring valley, ^vh, is ^vh '^0.1 pva ^ £ with 
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i the magnetic length. Vdonor is partly screened by the 
electron-electron interaction Vcc- The scppaied potential 
calculated in the Hartree approximationca'cS has an am- 
plitude approximately equal to the Landau level separa- 
tion, hujc- In such a potential, electrons (holes) populate 
in the 1th (A^th) Landau level, even in a filling factor 
of 2N {N = 1,2,---). 

In deriving the HDEQs, we take into account pro- 
cesses which are not included in the Hartree approxi- 
mation. We divide the electron-electron interaction into 
two terms: T4e = (Kse) + (Vee — (Vee))- The first term 
is the electron-electron interaction averaged with respect 
to the distribution of other electrons in the Hartree ap- 
proximation. (Vcc) and Vdonor give the aqremed poten- 
tial as calculated by the previous authors. Ej'ES' We intro- 
duce single-electron states in such a mean-field potential. 
The second term, on the other hand, gives scattering 
processes between the single-electron states. We con- 
sider such scattering processes only in the lowest order 

of Fec-(V^ee). 

2.1.1 Model of electronic states and processes 

Electronic states in the slowly- varying screened poten- 
tial consist of localized states with size ~ £vh around a 
potential hill or valley and extended states with energy at 
the center of the broadened Landau level. Each state has 
a level broadening due to the electron-electron scatter- 
ings. An estimate of the level broadening is F ^ 0.2hujc 
and that of the inelastic scattering length is ice ~ 
according to the estimation in §4.7 at Tc ~ Td (the def- 
inition of Tci is given in § p.l.2 ). Due to such electron- 
electron scatterings quantized energy levels of the lo- 
calized states are smeared out. We therefore label the 
localized states by a continuous momentum a which is 
related to a continuous energy in terms of the local elec- 
tric field. The same electron-electron scatterings destroy 
the phase coherence of extended states. In such situa- 
tion it is convenient to introduce wave packets with size 
^ £ee by superposing extended states. Here we divide 
the 2DES into regions so that each region contains one 
valley or one hill separated by equipotential lines of the 
fluctuating potential at the energy of extended states. 
We include each wave packet traveling along the bound- 
ary equipotential line in either the neighboring hill or 
valley region. All states are then labeled by a, i, and a 
where a — Na labels a Landau level with N = 1, 2, • • • 
the Landau index and a the spin, i labels a region, and 
a labels a momentum. 

Processes which transfer electrons between states con- 
sist of drift motions and scattering processes. Drift mo- 
tions perpendicular to the local electric field transfer 
electrons in extended states between neighboring regions 
(which are called d-processes) , giving the Hall current 
perpendicular to the macroscopic electric field. Scatter- 
ing processes we consider are due to electron-electron 
and electron-phonon interactions. In our model of scat- 
tering processes, a transition of an electron is restricted 
to be within the same level (diagonal in a) and no 
inter-Landau-level transitions are considered, since the 
transition rates are muciL. smaller than those of intra- 
Landau- level transitions .E2P In electron-phonon scatter- 



ings (which are called p-processes) we consider scatter- 
ings within the same region only (diagonal in z). In 
electron-electron scatterings, we consider the following 
two types. In the first type, one electron hops between 
neighboring regions, that is, between a hill and a val- 
ley, and another makes a transition within a region (h- 
processes). In the second type, both electrons make a 
transition within a region (s-processes). The temporal 
evolution of the occupation probability faia(t) of state 
aia is given by the Boltzmann-type equation: 



dfa 



dt 



dt 

^ foLia 

dt 



dfa 



dt 

^faia 

dt 



(2.1) 



2.1.2 Characteristic scales 

The lower critical current density jd and the lower 
critical electric field Ed are defined by the values at 
which the large-p^a; state becomes unstable. Values of 
the current density and the electric field are of the order 
of jci and Ed, respectively, in this paper. Similarly, val- 
ues of To are around Td, which is the value at jd in the 
large-pxx state. 

The length scales in the present system are the mag- 
netic length £, the average distance between a valley and 
a hill £vh, and £d defined by eEdici = hujc- £ci is shown 
in §5.3 to be the length scale of variations of Te and ^. 
The typical sizes are £ = 0.01/^m, ^vh ~ 0.1/J,m, and £ci ~ 
2fim at i? = 5T and Ed = 50V/cm, and £ < ^vh < £ci- 



The corresponding time scales are 



vh — £vh/^vh7 



and Td — £d/vd where v^h = eE^i^i /h is the group 
velocity at the local electric field between a valley and 
a hill i?vh = ^Wc/e£vh, and Vd = eEdi'^/h. The ratios 
between the time scales are Wc^vh = ^vh/^^ and 
^^cTci = £d/£'^ ^ iC. Finally the energy scale is htUc- 



2.2 Local Equilibrium 

We assume the local equilibrium in the scale of ivh 
and introduce the electron temperature Tc(ri, t) and the 
chemical potential fi{ri,t) of the ith region at the po- 
sition Vi ~ {xi,yi). In the local equilibrium faia{t) be- 
comes 



faiait) = fieaia - Ai(T-., , t) , To (r; , t)) , 



(2.2) 



where f{e,T) — l/[exp(e/A;BT) -I- 1] and Saia is the en- 
ergy of state aia. We neglect deviations of faia{t) due 
to the presence of electron and energy fluxes. With use 
of such averaged occupation probability, temporal evolu- 
tions of the system become irreversible. 

To be consistent with the local equilibrium in the scale 
of £vh , it is necessary that the inelastic scattering length 
£cc of extended states should be shorter than £vh- The 
estimation in §4.7 gives ice ^ ^vh when Tq ~ Td. Note 
that £cc becomes much longer than ^vh in the very-low- To 
region, in which the local equilibrium is violated. 

It is also necessary to have efficient electron and en- 
ergy transfers among different levels with different N and 
a in order to keep the local equilibrium. Therefore we 
must restrict our discussion to the case where temporal 
variations of macroscopic variables are slow enough to 
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keep the local equilibrium. 

2.3 Spatial Averaging 

The 2DES actually has a random configuration of po- 
tential hills and valleys instead of a periodic array. Due 
to such randomness, electron and energy densities, as 
well as electron and energy flux densities, fluctuate from 
place to place. We here take an average of these quan- 
tities over the area larger than the scale of fluctuating 
potential, £vh- Owing to this averaging, the equations 
become translationally invariant. To justify this aver- 
aging, it is necessary that £vh is much shorter than the 
length scale £var of variations of an d fi. This is approx- 
imately satisfied, since it is shown in § |3.3| that £var ~ ^ci- 
Using £vh ^ ^var, the difference of a quantity between 
neighboring regions is replaced with the differential in 

§1 

The averages of the quantities depend in general on 
the distribution function of the distance between a hill 
and a valley, which is mainly determined by the average 
£vh and the standard deviation. We neglect effects of 
the standard deviation in the calculation in the present 
paper. 

2.4 Other Assumptions 

We also make minor assumptions in order to simplify 
the HDEQs, although it may not be difhcult to go beyond 
them. 

(a) We neglect the change of electronic states (in the 
single-electron approximation) with the applied electric 
field since Ed ^ 0.1i?vh- 

(b) The screened potential depends on because the 
local electron density depends on To. We neglect such 
dependence by considering only the temperature range 
around Tci- 

(c) We neglect the broadening of energy spectra of 
extended states due to electron-electron scattering pro- 
cesses and the resulting decrease of the activation energy, 
which are of higher order in electron-electron interac- 
tions. 

(d) In addition to the valley-hill hopping processes 
which we have introduced above, elastic tunneling pro- 
cesses giving valley-valley and hill-hill transitions of an 
electron also contribute to the current along the electric 
field and consequently to pxx- We neglect such tunneling 
processes in this paper, for simplicity. 

(e) In calculating terms in the HDEQs, we assume that 
the change in energy of an electron at a scattering process 
Ae ^ji_,7r/cBTc (Ae eE^hi as shown in our previous 
papeirH') and that the lattice temperature Tl = 0, and 
neglect the energy dependence of the density of states 
and the spin splitting. 

§3. Hydrodynamic Equations 

3.1 Macroscopic Variables 

Macroscopic variables in the present system are the 
electron temperature Tc{r,t), the chemical potential 
IJ,{r,t), and the electric field E{r,t), with r = {x,y). 
We introduce a reference energy eni{r,t) defined by 
Em — i^Ni + eAf+ii)/2, where Sa is the energy of ex- 
tended states in the ath level, and the indices N and 



N + I are those of the two Landau levels in the both 
sides of fj,. The deviation of from Em is introduced and 
is denoted by fj,,- = ^ — Sm. The electric field is given by 
E = {Ex,Ey) = Vcm/e where V = {d/dx,d/dy) and 
e > 0. 

Tc and /ir are intensive variables corresponding to the 
energy density U and the electron density n, respectively, 
and are determined by the evolution equation of U and 
that of n, which are derived in the next two subsections. 
The increase of To produces the increase of pxx, and the 
change of /if at high To gives the deviation of pxy from 
the quantized value, as shown in eq.(4.17). 

3.2 Equation of the Charge Conservation 

The electron density averaged over regions within area 
S satisfying £vh < S^/"^ < 4i is 



(3.1) 



Here the integration over a continuous momentum a is 
denoted by the summation, for simplicity. The drift and 
hopping processes give electron transfers between neigh- 
boring regions. Since the change of n is given only by 
electron transfers between the area 5* and the outside 
through the boundary L, we have 



dt 



(3.2) 



with Ji the number of electrons going from the ith re- 
gion to the outside of S per unit time. Introducing the 
averaged electron flux density _7„, we have 



ndL 



V-i„d5 , 



(3.3) 



where n is the unit vector perpendicular to the bound- 
ary L directed to the outside, and the two-dimensional 
Gauss theorem is used in the second equality. Since the 
variation of V • j„ within S is negligible, we obtain the 
equation of the charge conservation: 

^^_-V.,„. ,3.4, 

The electron flux density j„ is given by 

where j'jj and j'Jj are the electron flux densities due to 
drift and hopping processes, respectively. 

3. 3 Equation of the Energy Conservation 

The energy density due to the thermal activation, aver- 
aged over regions within area S satisfying ^vh ^ S^^^ 
4i, is 

U{r,t) = i ^ (£„,, - /i,)(/a.a - /°.J , (3.6) 
a.i,a 

where fii = //(/■;, t) and /"j^ is the occupation probability 
in the local equilibrium at Te = 0: f°^^ = f{£aia - Pi,0)- 
Terms with Saia > fJ-i express electron excitations, and 
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those with Saia < fJ'i express hole excitations: 



/ ocia 
otia f^i ) /as 



fai 



< Mi 



U is considered as a function of the electron dis- 
tribution described by and faia- In dU/dt, terms 
with d/j-r/dt are neglected since Hr changes much slower 
in time than faia- Then dU/dt is given, using the 
Boltzmann-type equation eq.(2.1), as 



dU_ 

'at 



(du\ 



'dt 



dU_\ fdU\ /dU 



at y, 



hi \ / h2 

Here we have decomposed the term due to h-proccsscs 
into two terms: hi and h2 denote the contribution 
from hopping electrons and that from electrons scattered 
within a region, respectively. 

The drift and hopping processes give energy transfers 
between regions. The transferred energy \eaia — Mil is 
£ra " Mr electron excitations and /ir — Era for hole 

excitations with Era = £a — £m- Although Saia — £m of 

the initial and final states in the hopping processes has 

fluctuations around Era due to the fluctuating potential, 
such fluctuations are canceled by those in the energy of 
the scattered electron in the calculation of the h2 term 
(see below). The s-processes also give energy transfers 
between regions. The change of U due to such energy 
transfers is given by 



dt 



dU 
dt 



dU 
'dt 



d \ / hi 

following the procedure in the previous subsection. Here 
the averaged energy flux density jjj is given by 



Ju 



Ju 



(3.10) 



where ju, ju, and are components due to drift, hop- 
ping, and scattering (without hopping) processes, respec- 
tively. 

The 2DES gains an energy from the electric field 
through the energy change As of the scattered electron 
in the h-processes. Neglecting the fluctuations in As (see 
above), Ae = {—e)E-Ar with Ar the change in position 
of the hopping electron. By summing over all hopping 
processes within S, we obtain 



/dir 
yat 



) ={-e)E.jl 

/ h2 



(3.11) 



Electron-phonon scatterings give changes of energy 
within each region and result in the energy loss after 
the statistical averaging when the lattice temperature is 
lower than T^: 



dU 
'dt 



= -Pt, 



(3.12) 



where Pl > is the energy loss per unit area per unit 
time. 



Then we obtain 

dU 
'dt 



V -ju + {-e)E-j^-P, 



L , 



(3.13) 



(3.7) where we have used E ■ j'^ = 0. 



3.^ Boundary Conditions 

We impose boundary conditions at the edge of the 
2DES in solving the HDEQs. In the case of the boundary 
between the 2DES and the depleted region, they are 



j„ • nb = , i[/ • nb = , 



(3.14) 



where n\, is the unit vector perpendicular to the bound- 
ary. Another boundary condition is that the total cur- 
rent is given by the external circuit. 

§4. Evaluation of Terms in the HDEQs 

4.1 Dimensionless HDEQs 

Here we introduce dimensionless units in which the 
length and the time are scaled by £ and ui~^, respectively. 
Dimensionless variables are defined as f = r/£, V = ^V, 
i = ojct, and E = eEE/htOc- We use TiUc to scale energies 
and introduce = k^T^/TiuJc, Mr = Mr/^i^o and ^aia = 

£aia/'fVjJc- 

In the present units the equation of the charge conser- 
vation becomes 



dh - - 



(4.1) 



with h = i'^n and j„ = {£/u}c)jn- The equation of the 
energy conservation becomes 



^ = -t-ju-E.j^-Pr. , (4.2) 

with U = {eynio^p, ju = {£/hiol)ju, and Pl = 
{e/h^l)Pi.. 

4-2 Time Derivatives of Electron and Energy Densities 
The electron density is given in the local equilibrium 

by 

pOO 

n{r, t) = 2 dsp{e)f{s - m, T,) , (4.3) 

Jet, 

where p(e) is the density of states per spin per unit area 
and Eb is the energy of the lowest state. Using the as- 
sumption that p{e) is constant: p = l/{2Tr£'^hujc), we 
obtain n(r , t) = 2p{p — e\y) = 2p{pr + Em — Eb) as long as 
/(sb — M)^e) ~ 1- In the dimensionless units, we have 

dt TT dt 

Similarly, the energy density is given in the local equi- 
librium by 

/oo 
d£p{£){£-p) 
-00 

x[/(E-M,Te)-/(E-M,0)] . (4.5) 
Since p is assumed to be constant, U{r,t) does not de- 



pend on p.^, and we have 

dU _ dU_m 
'dt ~ dT^~dt 

with the dimcnsionless specific heat Cc = ■^T,^ 
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extended states and is given by 



(4.6) 



a.i.a.b 



{iaia - £aib)faiai^ - fa^b)W^^^_^^f, , (4.7) 



4-3 Energy Loss 

The energy loss due to electron-phonon scatterings is 
given by 

LX.L 

where S = S/£^ and W^^^^^^i^ is the number during 
the time uj^^ of electron-phonon scatterings, in which an 
electron makes a transition aia — > aib, and a phonon is 
created or destructed depending on the sign of Saia—Saib- 
We evaluate Pl when the lattice temperature Tl — 0. 
Here we assume that \eaia — £aib\ ^ T^kBTg, which is 
approximately satisfied since \saia — £aib\ ~ eE^iiiy 
O.lhujc and irkBTe ~ TllUc at To = Td from eq.(5.3). 
Using this assumption, we have /aw(l - fmb) « /aw(l - 
faia) = -kBTcdfisaia ~ fJ.,0)/ds, and then obtain 



Pl = C„T, 



(4.8) 



with 



a.i.a 



(-l)^(ea. 



M,0) 



it 



(4.9) 



We here neglect the weak aia dependence of W^^^_^-f^ 
(Wma^ib depends more strongly on Satb - £aia through 
the overlap integral between the states). Using also the 
assumption of no energy dependence of the density of 
states, Cp becomes independent of Tc and /ir, and is 
given by 

Cp = i^(£„., -£„,)<„ . (4.10) 



4-4 Drift Components of Electron and Energy Fluxes 

The electron flux density, averaged in a macroscopic 
scale, due to drift processes, is given by 



~d 



7ia "> 



(4.11) 



-d 



where is the averaged electron flux density carried 
by electrons in the level a. The energy flux density, on 
the other hand, is given by 



3u 



i(ea>M) 



a(eQ<M) 



(Mr 



l^r)3na 



3 na } 



(4.12) 



where j'^^i is the value of j*^^ when the level a is com- 
pletely filled by electrons. 

Since j^^ is an averaged flux density, it is only from 



~.d 

J na 



faJnal 



(4.13) 
(4.14) 



fa fi^ra /^r,Tc) . 

It was shown by MacDonald et that the local elec- 

tron flux density of the filled Landau level due to the drift 
motion in the local electric field Eioc is independent of 
the Landau index (therefore denoted as j^l"'^) when E\oc 
varies in the scale » £, and that j^'l"'' — {e/h)E\oc x e^, 
with Bz the unit vector along z {B — Be^ with B > {)). 

In the dimcnsionless units, j^^ = ^Eioc x e^- By 
averaging this equation over 5*, we obtain 



~.d 



X 



^{Ey, -E^) , 



(4.15) 



and 

given by 



~.d . ~d 

jj^i- The electron flux density j„ is then 



~.d 

3 n 



^xy{-^y: Px 



fa 



(4.16) 
(4.17) 



0, a^y = with N the number of filled 



When fir 
Landau levels. 

The energy flux density, on the other hand, is given 

by 



~d 

3u 



^xy{,Pyi Px) : 
a{ea>ti) 



2ir / y 

q(£q<Ai) 



{fir " £ra)(l " fa) 



(4.18) 



(4.19) 



When we consider only the two Landau levels near /x, as- 
sume fir — 0, and use the assumption of no spin splitting, 
we obtain u^y = -^fi with 

-1 



A 



exp 



1 



(4.20) 



4-5 Hopping Components of Plectron and Pnergy 
Pluxes 

The electron and energy flux densities, averaged in a 
macroscopic scale, due to electron-electron scatterings 
(h-processes) are given by 



~.h _ 1 

~h _ 1 
3v--g 



4d > 



(4.21) 



a^f3^i,k,l,a,b,c,d 



J2 - Ar)(r/ - H) J'^i:^ , (4.22) 



where S = S/^ and 

JpkiPdd = faia{l-faib)ff3kc{l-ft}ld)Wpl''^Zld ' (4.23) 

with Wpl'^'^l^ the number of scatterings during the pe- 
riod Lo^^ ■ Regions k and I are the nearest neighbors. 
For hole excitations {erp — /2r < 0), the above formula of 

3ij is consistent with the picture that each hole excita- 
tion carries the energy of fir — irp, since the hole flux is 
opposite in sign to the electron flux. 
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~h 



We rewrite j„ as 
1 

2^ 



a.J3,i,k,l,a,b,c,d 



jaib — na \ 



(4.24) 

and Ju similarly, where changes of labels are made in 
the term with JpijZ*,]^^ '■ ci ^ b and kc ^ Id. Since the 
two terms are forward and backward transitions among 
the same set of states, Jplc^l^ — J/^mZIZ ^^'^ "^^^ 
number (during cj^^) of electrons which hop from c in fc 
to d in /, and is induced by V/i and VTo. In the first 



order of V^i and VTe, JptcZld 



j: 



aib — ^ia 



jata- 



*ld 



jatb — >ia 
•Jpid-^kc 



-T~ 



is given by 



jaia — nb 



ll jaia — 'ib 
\ '^f3kc^ld 



(4.25) 

where jlk = p.{rk,t) /fiujc, fk = fc{ru,t), and Jp^Ztd 
in the right should be evaluated in the absence of 
V/x and VTc. Higher-order terms are negligible when 
\Vii%h/kBTc < 1 and | Vre|4h/2Tcfe < 1. These 
conditions are approximately satisfied at £^ = E^i and 
To = Tci- since the estimated value of E^i in eq.(5.6) 
and that of Tci in eq.( |3.3[ ) give | V/zl^vh/^BTc ~ 0.1 and 
I VTc|^vh/2T'oTc ^ 0.1. We also have made an approxi- 
mation that Spkc - l-ik^Sfild - fM ~ £i-i3 - Mr- 
We then have 



3n 



— y {ri 

a,f3,k,l 



f-k) 



X {fi - f-k) ■ [-f-'Vil + (e> - Ar)VT-i] 



E 

i^a^b,c,d 



jaia — >ib 
•-'pkc^ld 



(4.26) 



Here we make the integration over 9ik, the polar coordi- 
nate of ri - rk = rikexp{i6ik) = (xjk.yik)- The depen- 
dence on ri — rk appears through Jplc^id ^^'^ ~ ''^k, 
but the 0ik dependence appears only through f; — f fc, 
since 'Et,a,b,c,d JplcZld has no dependence on 0ik. We 
assume the uniform distribution of 9ik, in which we have 
{xikVik) = and ((xzfc)^) = {(yikf) = {rikf /I, with 
(• • ■) = ^ /o ' ' ' d^;fc- We also make an approximation 
that rik « ^vh- h h 

Then we obtain the formulas of j!„ and jjj: 

~h 

3n = 



3u 



L 



^(-f-iV/i) +L„TVr-i , (4.27) 
(~f-^Vji) + , (4.28) 



with 



Ln^i = Lq , LnT = Lij^i = Ll , LuT = L2 , (4.29) 



i: 



T vh \ ^ /- vm Jaia^ 

' ~ ^ — Mr J ■Jf3kc^ 



ib 

Id ) 



a,j3,i,k,l,a,b,c,d 



(m = 0,1,2) , (4.30) 



where V/i = eE + V/ir ~ e.E and in dimensionless units 
V/i w E. 



Here we assume \eaia — £aib\ ■C tt/cbTo. Then we 

have /aia(l ~ faib) ~ faia{^ ~ jaia) = ~kBTcdf{£aia ~ 

/i,To)/5£ and f/skcfpid « ffS- We neglect weak depen- 
dences of W^fcc^/d on a,i,a, and /?, for simphcity. We 
then obtain, by performing the summation over a, i, a, 



iafeE(£~r/3-Ar)"//3(l~/;3) 



/3 



Here 



vh \ Ti^aia— * 

<-h = — 2^ Wpkc^ 



(4.31) 



(4.32) 



kJM.cA 



with ^vh = ^vh/^j has no dependence on To and /ir. 

When we consider only the two Landau levels near /i, 
assume /ir = 0, and use the assumption of no spin split- 
ting, we obtain L„^ = Chre/i(l - /i), T„t = Luii = 0, 
and ic/T = \Ch_Tcfi{l — fi). In the present case L„t = 
because energy levels of extended states are symmetric 
with respect to /i and electron fluxes induced by VTo are 
canceled between the upper and lower Landau levels, and 
Lu^J. = because electrons and holes move in opposite 
directions by hopping processes in the presence of E, in 
contrast with drift processes, and energy fluxes due to 
such hopping processes are canceled between the upper 
and lower Landau levels. The electron and energy flux 
densities in this case are given in terms of the electrical 
conductivity axx and the thermal conductivity Kh as 



3n = -(^xxE , 



3u 



-KhVTc , 



(4.33) 
(4.34) 



with 



a.x/(eV?i) = Ch/i(l-/i) , (4.35) 



(4.36) 



4-6 S-process Component of Energy Flux 

The energy flux density, averaged in a macroscopic 
scale, due to electron-electron scatterings without hop- 
pings, is 



3u 



where 



E i^aib — £aia){'f'i — fk) Jf3kc~*^'' 



2S 



a.,l3,i,k,a,b,c,t 



(4.37) 



Taia—>-ib \ /i £ \Tjra.ia — >ib 

■Jpkc~>kd - Jaia\^ ~ Jaib)Il3kc\^ - jpkdjyy pkc^kd ■ 

(4.38) 

The factor 1/2 in the formula of corrects the double 
counting: i = 1, fc = 2 and z = 2, fc = 1, for example. 
Following the derivation of j'^j in the previous subsection, 
we obtain 



(4.39) 

when I VTc I ^vh /To ^ 1, with the thermal conductivity 

'*s — ^^2^2 E Jpkc^kdi^aib—^aiaYifi—fkY- 

c a.p,i,k,a,b.c.d 

(4.40) 

Here we assume that \eaia~£aib\ ^ TrfceTo and neglect 
the weak a, (3, a, c dependence of 1^^^"^^^ ■ We restrict 
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the calculation to the case of /ir = and no spin splitting. To at Ed in the high branch is Td, and is given by 
Then kg at ^i- = 0, which is denoted by KsO, is given by 



e 



i.k.h.d 



kso has no dependence on 



(4.41) 



4-.7 Inelastic scattering time 

Here we calculate the inelastic scattering time Taia of 
state aia due to electron-electron scatterings. Consider- 
ing s-processes, Tata is given by 

= X! ~ fatb)f/3kc{l - f fUkdW^l^^Ztd ■ 

'^'^^"^'^ fS,k.b,c,d 

(4.42) 

Following the evaluation of ks, we assume \£aia — £aib\ ^ 
TrfesTe and neglect the /9, c dependence of Wj^^'^Zlkd ■ 
Then we obtain, when /ii- = 0, 

1 



^ib 
kd 



(4.43) 



k,b,d 



An estimation gives 1/iOcTaia ^ 0.2 for an energy £aia — 
Em ~ huJc/2 at /ir = 0, B =5T, and = Td (the value 
of Td is given in eq.( ^.3|) ). Then the level broadening F 
is estimated to be F = h/Taia ^ 0.2?iWc. The inelastic 
scattering length defined by £cc = ''aiai'vh is estimated 

to be ico/( = {LUcTaia)eE^i,i/huJc 1- 

§5. Estimation of CoefRcients and Scaled HD- 
EQs 

5.1 Uniform Steady States at Even-Integer Filling Fac- 
tors 

We here consider the simplest case of uniform steady 
states at even-integer global filling factors and obtain 
the formula of the critical electric field Ed as well as 
the value of Td- In the present case we have fi^ — 
everywhere in the sample. We consider only the two 
Landau levels near /i, neglect the spin splitting, and then 
use the simplified formula of (Jxx = Ch/i(l ~ /i)- We 
will give only a summary of tho-ccsults, which we have 
obtained in our previous paperBiJ) The equation of the 
energy conservation in this case gives the balance of the 
energy gain Pq and the energy loss Pl: 



{Po-) 



rE- 



= CpT, ( = Pl 



(5.1) 



Pq and Pl are plotted as a function of To in Fig. 1(a). 
Points of intersection of the two curves, PG(Pe) and 
Pl{Tc), give Te in stationary states. The number of 
points of intersection increases from one to three as in- 
creasing E through a critical value. This critical value is 
Ed, and is given by 



Ed = (Cp7/Ch)i/2 (7 = 2.2334) 



(5.2) 



A bistability appears above Ed (one of three in the mid- 
dle corresponds to an unstable state) . The present bista- 
bility originates from the existence of the activation en- 
ergy hujc/2 of Pq. Tc and axx in stationary states as a 
function of E/Ed are plotted in Fig. 1(b). The value of 



Td = 0.324 . 



(5.3) 



Fig. 1. (a) The energy gain Pq = axxE^ at the critical electric 
field Eel and the loss Pl = CpTe (both divided by Cp) as a 
function of the electron temperature Te. Points of intersection 
give Tc in stationary states, (b) Te and axx/Ci^ in stationary 
states as a function of E/Eci- 



5.2 Estimation 0/ Coefficients 

Rough estimates of the coefficients in the KDEDs at 
B = 5T have been given in our previous papersE2HiH' and 
are summarized here: 



Cp = 1.0 X 10" 



Ch = 2.5 , = 0.016 



(5.4) 



The ratio of the two components of the thermal conduc- 
tivities, Kh and Kg, is then estimated to be 



(5.5) 



when /ir = and Te = Td- The thermal conductivity 
due to hopping processes is much larger than that due 
to electron-electron scatterings without hoppings. In the 
following we neglect ks- 

The above estimates of Cp and Ch give theoretical es- 
timates of Ed and — jC^, the saturation value of 
axx at large To (only the two Landau levels near fi are 
considjeied). The theoretical estimates and the observed 
valuesErEJ' are as follows. 



Theory(P = 5.0T) : E'd = 3 x 10^^ 
Rcf. |(P = 4.7T) : Ed=9x 10"^ 
Ref. Mb = 3.8T) : Ed^S,^ 10"^ 



xx~ 



0.6 , (5.6) 
0.02 , (5.7) 
0.02 . (5.8) 



The discrepancies between the theory and the experi- 
ments will be within the limitation of accuracy in the 
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theory, since a large ambiguity in the estimation of i 
produces large errors in the estimation of Cp and Ch- It 
may be convenient to consider Cp and Ch as adjustable 
parameters when spatio-temporal evolutions are com- 
pared between the theory and the experiment. 

5.3 Scaled HDEQs 

The energy gain and the energy loss in the HDEQs 
are very small in magnitude in the units we have used: 



Pi 



10-5 and Pl ~ CpTci - 10" 



m 
V 



2.1.2 



according 

to the estimation in the previous subsection. Therefore 
we switch the units to the ones in which the length and 
the time are scaled by id and Td, respectively, defined 
Scaled variables are defined as r = r/ld, 
iV, i = t/rci, and E = eiciE /TiiOc- Ratios be- 
tween the present scale and the previous scale are ^ci/^ = 
{Eci)~^ and TdWc = {Eci)~^- Since E — E/Ed, the en- 
ergy gain and loss become of the order of unity in the 
present units. 

When fi,- — and only the two Landau levels near fj, 
are considered for hopping processes, the HDEQs in the 
new units become 

1 dfir 



(5.9) 



- a,,,E^ - C^j-^f, , (5.10) 



with 



On = (4l/Wc)j„ = -^xxE - G^y[E X 6;,) , (5.11) 

3u = {£d/huj^c)3u = -KhVTo + u^y{E x e,) . (5.12) 

The scaled HDEQs show that the typical length and time 
scales of variations of and are td and Td, respec- 
tively. 

§6. Current Distribution at Even-Integer Filling 
Factors 

In this section we illustrate how to apply the HDEQs 
to the calculation of spatial variations of Tq. We consider, 
as an example, the current distribution in a steady state 
in the higher- To branch in the middle of a long sample, in 
which spatial variations occur only perpendicular to the 
current. We restrict our calculation to the simple case 
where the global filling factor is 2N (N = 1, 2, ■ • •) and 
To is low enough that only the two Landau levels near fi 
contribute to jjj, j^, and jfj. 

6. 1 Electrostatics and Estimation of fi-^ 

In this subsection we derive the equation relating E 
and /ij.- We show here that //j. ^ {EdY 10"^. 

We define the two-dimensional local charge density by 
'y\oc{f) — —e{n\oc — «o) with rtioc the local electron den- 
sity and rto the uniform electron density corresponding 
to the global filling factor 27V. The local charge density 
is divided into two components: crioc = cr -\- with a 
long-ranged variations due to the Hall polarization and 
(Jsr short-ranged fluctuations induced by density fluctu- 
ations of ionized donors. The local electric field E\oc is 
also divided into two components: £^ioc = E + Esr with 
E the macroscopic electric field consisting of the Hall 



field induced by a and the uniform field along x and £^sr 
the short-ranged fluctuations due to ionized donors and 
(Tsr- The long-ranged variations given by a and E are 
assumed to occur only along y {0 < y < W). 

The dimensionless charge density a — aP' je is given 

by 



7 Mr 



The first term of a was given by MacDonald et alB 
It is the charge density in the presence of spatial mod- 
ulations of the electric field when the lower iV Landau 
levels are fully occupied and the other higher Landau 
levels are empty, and it comes from the polarization of 
wave functions of the occupied states. The second term 
is the charge density due to the occupation of electrons 
or holes, which is described by /Zr when averaged in a 
macroscopic scale. 

The charge density a is also given by the Gauss law 
relating a and Ey: 

e'/ee f^/'^^, 2a{y')sgn{y-y') 
Ey{y) = — / dy — , , (6.2 

where e is the dielectric constant and £zc is the cutoff 
length corresponding to the finite thickness of the elec- 
tron density perpendicular to the layer. Here sgn(?;) ~ 1 
when y > and sgn(?/) = — 1 when y < 0. 

Now we estimate the order of magnitude of fir- From 
eq.(^), we have 



dEy 

dy 



huJr 



4CT(y) 



W/( 



dy' 



2aiy')\y-y'\ 



[{y 



y'Y -f z2]3/2 



huJc Zr 



(6.3) 



where we have assumed that does not change signif- 
icantly in the scale of z^. Then we have (T(j/) ~ dEy/dy 
since e^/ef ^ fiLUc and z^ ^ \. Combining this with 
eq.(O), we obtain jl^ ^ dEy/dy ^ Edited = (£^ci)^ ~ 
10"^ 

6.2 Macroscopic Variables and Equations 

It was shown in the previous subsection that jl^ ^ 
[Ed)^ ^ 10~5 when the system is uniform along x and 
the global filling factor is an even integer. Since /ir ^ 1 
(/ir ^ 'hijJc) and jl^ <^ (/ir k^Tc) , we neglect 
jlr in the HDEQs and use eqs.(5.9), (5.10), (5.11), and 
(5.12). Macroscopic variables Tc{y) and Ey{y) are ob- 
tained by solving the two HDEQs, and fJ.r{y), on the 
other hand, is determined by eqs.( |6.1| ) and ( |6.2| ) in the 
previous subsection, using the obtained Ey{y). In the 
present case, E^ does not depend on y, either, since 
dEx/dy = dEy/dx = 0. The constant value of E^ is 
determined by the total current through the system. 

6.3 Equations for Tc{y) and Ey{y) 

The equation of the charge conservation in the present 
case is djny/dy = 0, and the boundary condition is jny = 
at sample edges {y — 0, W). Then we obtain 



Jr. 



-{Pxx 



Er, 



(6.4) 
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jny — —O-yyEy — (JyxEx — , (6-5) 

for < y < T4^, with = axx/{<^lx + '^xy)- The equa- 
tion of the energy conservation in the present case is 

^R^{f,)^f, = -nit)-^f, + Fit) , (6.6) 

Vit) ^ ^Exfiil - fi)f-\ (6.7) 

F{f,) = ~p-^El + Ch7"'Te . (6.8) 

Here we have used eq.(6.5) to ehminate Ey. The bound- 
ary condition is juy — a.t y ~ 0, W, and is given by 

^h^Te - -yiEx . (6.9) 
oy 

6.4 Current Density Modulation 

The above boundary condition eq. (|6.9[) shows that To 
depends on y due to a nonzero Ex- This modulation 
of Tc gives a modulation of jx = {—e)jnx, since pxx 
in eq.(6.4) depends on Tc (when jjxx is a constant, on 
the otheiLJ|iand, jx becomes uniform as pointed out by 
ThoulessE3). The present current density modulation 
originates from the drift component of the energy flux 
perpendicular to E (the second term in eq.(5.12)). This 
energy flux has nonzero y component when E^. ^ and 
produces the modulation of To in eq.( |6.9| ), and conse- 
quently the modulation of jx ■ 

Since the equation for is nonlinear, numerical cal- 
culations are necessary to obtain the solution, which will 
be performed in the future work. Here we instead discuss 
the qualitative profile of Tc and jx by regarding eq.(6.6) 
as an equation of motion of a particle with the "position" 
To as a function of the "time" y. In this mechanical ana- 
logue, Kh is the "mass" , 77 is the "friction coefficient" , and 
F is the "force" . The corresponding "potential" around 
the higher- To stationary point is plotted schematically 
in Fig. 2(a). The boundary condition eq. (|6.9| ) gives the 
"velocity" of the particle at the "time" y = 0, W/£ci- 

When W is large, the particle spends the most of its 
"time" at the top of the "potential". This means that 
To(j/) and jx{y) are approximately uniform in the middle 
of the sample, as shown schematically in Fig. 2(b). The 
length scale for the relaxation of Tc{y) and jxiy) from 
their deviations at the edges is of the order of ^ci 2/im. 

In the QHE regime with negligible pxx , MacDonald et 
alB^ have shown that jx{y) is enhanced near the both 
edges, because of the edge electrostatics in the 2DES. 
When Pxx is large, on the other hand, we find in Fig. 2(b) 
that jx{y) in one edge is enhanced, while jxiu) in the 
other edge is suppressed. This comes again from the 
boundary condition eq.( |6.9[ ) stating that the derivative 
of Tc{y) has the same sign in the two edges. 

§7. Conclusions and Discussions 

We have derived hydrodynamic equations (HDEQs) 
for quantum Hall systems in the regime of large energy 
dissipation. We have shown that the strong magnetic 
field in the present systems produces a component of the 
energy flux which is perpendicular to the electric field. 
The presence of this energy flux is the most important 
feature of our HDEQs for the quantum Hall systems in 



Fig. 2. (a) "Potential" in a mechanical analogue as a function of 
the "position" Tc in the vicinity of the higher- Te stationary point, 
(b) Schematic spatial variations of the electron temperature Tc 
and the current density jx in the sample (shaded region) when 
jx > 0. 



the high-electron-temperature regime. This energy flux 
is absent in the equation proposed by the previous au- 
thorsO The HDEQs with this unique energy flux have 
predicted a new type of current density modulation. 

Another example of the application of the HDEQs is 
to investigate spatial evolutions of To in the current di- 
rection (along x) when a spatial transition occurs be- 
tween the small pxxSt&te and the large pxxState along x. 
The observed spatial evolutions of pxx in such transitions 
have exhibited a divergence of the r elaxat ion-iongt h as 
the current density approaches a critical value. tj'lliHl3) In 
our previous paper we have applied the previous one- 
dimensional HDEQ to this problem and shown that the 
calculated relaxation length also exhibits a divergence 
qualitatively similar to the-observed one. 

A hydrodynamic modelES has been propos£|d_to ex- 
plain the observed longitudinal voltage steps.EIlES' It is 
argued in this model that electron-hole pairs created at 
inter-Landau-level scatterings near a charged impurity 
have common features with vortices created when a clas- 
sical fluid flows past an obstacle. The hydrodynamic 
equations in the present paper, on the other hand, de- 
scribe macroscopic spatio-temporal variations, and their 
typical length scale is much larger than the size of a sin- 
gle electron- hole pair. 
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